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Abstract 

We consider the model of strongly-correlated system of electrons described by an extended 
Falicov-Kimball Hamiltonian where the stability of some axial and diagonal striped phases was 
proved. Introducing a next-nearest-neighbor hopping, small enough not to destroy the striped 
structure, we examine rigorously how the presence of the next-nearest-neighbor hopping anisotropy 
reduces the 7r/2-rotation degeneracy of the diagonal-striped phase. The effect appears to be similar 
to that in the case of anisotropy of the nearest-neighbor hopping: the stripes are oriented in the 
direction of the weaker next-nearest-neighbor hopping. 

1 Introduction 
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{ Experimentally, the presence of the striped structures in materials exhibiting high-temperature super- 
conductivity is well-known pfl E] • Theoretical investigations of the stability and properties of striped 
^ phases are performed mainly with approximate methods. Usually, the Hubbard-like or t-J-like models, 
^ >as the standard models applied to explain the effects in strongly correlated systems, are considered 
£h for this purpose (see for example a review [3]). Unfortunately, due to the tiny energy differences be- 
tween the energies of compared phases, the methods that are applied to study the stability are not 
reliable. That is why some simpler models (see Refs. 01130111) are considered in the investigations of 
^ striped-structures stability and properties. 

In the recent paper jS] we addressed the issue of striped-order formation in the systems of strongly- 
•»H .correlated quantum particles described by extended Falicov-Kimball models. We showed that at half- 
rN [filling and in the strong-coupling regime some axial and diagonal striped phases are stable. In corn- 
ed jparison to the standard spinless Falicov-Kimball model (like in Ref. [3 El) the model in Ref. [S] was 
augmented by a direct Ising-like nearest-neighbor (n.n.) attractive interaction between the immobile 
particles, in order to allow for segregated phases. Changing the intensity of the Ising-like n.n. interac- 
tion, the system is driven from a crystalline (chessboard) phase to a segregated phase, via quasi-one- 
dimensional striped phases. This conclusion was reached for two kinds of hopping particles: fermions 
and hard-core bosons. 

The obtained results enable us to continue the study of the properties of striped phases. One of the 
interesting questions, the influence of n.n. anisotropy of hopping intensity on axial striped phases, was 
investigated by means of the Hartree-Fock method in [HIE]- in Ref- at the regime where stripes 
are stable, we have proved rigorously that for both systems, of hopping fermions and hard-core bosons, 
an arbitrarily small anisotropy of n.n. hopping orients the axial striped phases in the direction of a 
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Figure 1: The phase diagram obtained in [H], containing a diagonal-striped phase Sdd, at /i = 0, in 
the plane {uj,e). In this paper uj is given by W = —It 2 + ut 4 , see Q, and e = 0. The representative 
configurations (up to translations) of phase Sdd are shown on the right. For more details see [S]. 

weaker hopping. We noted also the tendency of the phase diagrams for different statistics to become 
similar, even for a weak anisotropy of n.n. hopping. 

The analogous, arising naturally, question is how the anisotropy of the next-nearest-neighbor (n.n.n.) 
hopping influences the degeneracy of diagonal-striped phases. To answer this question, in this paper, 
we use the same techniques as in j!2j . Specifically, we investigate the influence of n.n.n. hopping on 
the phase Sdd, whose stability was proved for fermions in JHj (in Fig. [T]we reproduce the phase diagram, 
and show representative configurations of the phase Sdd)- According to the state of art, that is possible 
only in the strong- coupling regime and at the half-filling. 

In the next section we present the model and discuss some of its basic properties. In Section 3, we 
introduce a strong-coupling expansion of the ground state energy (effective Hamiltonian) . In Section 
4, we carry out the analysis of phase diagram due to the truncated effective Hamiltonian. Finally, we 
make conclusions in the Summary. 



2 The model and its basic properties 

The model we consider here, is an extension of the spinless Falicov-Kimball model ^Uj; which in turn 
is a simplified version of the one introduced in [13 . The system is augmented by a direct Ising-like 
interaction between immobile particles. Hopping particles are allowed to hop not only between n.n. 
sites (as is considered usually), but between n.n.n. sites as well. Moreover, the n.n.n. hopping intensity 
depends on direction. The Hamiltonian of the system is of the form: 

H = H FK + V, (1) 

t + ^2 ( c t c y + C y °x) — t- ^2 { c t c y + C y °x) 

(x,y)2,+ (x,y)2,- 

(2) 



h fk = -t 22 ( c t c v + 4 c *) - 

1' 



2 



(3) 



The underlying lattice A consists of |A| sites denoted x, y, . . ., and having the shape of a a/|A[ x a/|A| 
torus. The sum Yl{ x ,y)i stands for summation over all n.n. pairs, while the sums ^2/ XiV \ 2 , ^2{ xy ) 2 _ 
for summation over all n.n.n. pairs oriented in the direction with slope +1 (+) or in the direction with 
slope —1 (— ), respectively. 

The subsystem of quantum hopping particles (here-after called the electrons) is described in terms of 
creation and annihilation operators of an electron at site x, c x , c x , respectively, satisfying the canonical 
anticommutation relations. The total electron number, N e , is the eigenvalue of Y2 x c t c x, an d the 
corresponding electron density is p e = N e /\A\. Although the electrons do not interact directly with 
each other, there is an on site interaction with the localized particles that brings energy U if a site is 
occupied by the both kinds of particles. 

Since the site-occupation-number operators of classical immobile particles (here-after called the ions) 
commute with Hamiltonian ((TJ), the subsystem of ions can be described by a set of pseudo-spins {s x } x£ a 
(the ion configurations), with s x = ±1 (1 for the occupied site and —1 for the empty site). The total 
number of ions is Ni = J2 x (s x + l)/2 and the ion density is pi = A^/|A|. In our model the ions interact 
directly via a n.n. Ising-like interaction V, which is isotropic. 

The Hamiltonian Hfk is the well-known spinless Hamiltonian of the Falicov-Kimball model. A 
review of rigorous results and an extensive list of references concerning this model can be found in 

We shall study the ground-state phase diagram of this system using the grand-canonical formalism, 
i.e. we consider 

H (ii e , p^ = H o - p e N e - piNi, (4) 

where p e , pi are the chemical potentials of electrons and ions, respectively. Let E s (p e , p{) be the ground- 
state energy of H(p e ,Pi), for a given configuration S of the ions. Then, the ground-state energy of 
H(p e ,pi), E G (p e ,pi), is defined as Ec(p e ,Pi) = min {Es (p e , Pi) '■ S}. The minimum is attained at 
the set G of the ground-state configurations of ions. 

Applying unitary transformations to H (p e , p^) we can restrict the range of the energy parameters. 
Firstly, using the hole-particle transformation for ions, s x — > —s x , one finds that the case of attraction 
(U > 0) and that of repulsion (U < 0) are related by this transformation: if S = {s x } x< z\ is a ground- 
state configuration at (£,£+,£_, p e , Pi,U), then — S = {— s x } xe \ is the ground-state configuration at 
(t,t+,t-, p e , —pi, —U). Consequently, one can fix the sign of the coupling constant U without any loss 
of generality. We choose U > 0. We express also all the other parameters of the Hamiltonian (jjj) in the 
units of U, i.e. formally we set U = 1, preserving previous notations. 

Secondly, applying the hole-particle transformation for electrons, i.e. c x — > e x c x (where e x is equal 
to +1 on the even sublattice, and —1 on the odd one), and the hole-particle transformation for ions, 
one finds that if S is the ground-state configuration at {t,t + ,t_, p e , pi), then — S is the ground-state 
configuration at (t, —t + , —t-, —p e , ~Pi)- Thus, it is enough to consider only one sign of n.n.n. hopping. 
We shall consider the case of positive n.n.n. hopping intensities: t+,t- > 0. 

Thirdly, using another hole-particle transformation for electrons, c x — > c x , and the hole-particle 
transformation for ions, we obtain that if S is the ground-state configuration at (i, p e , p^), then 

—S is the ground-state configuration at (— t, —t + , —t_, —p e , ~Pi)- Applying consecutively the two joint 
(with respect to electrons and ions) hole-particle transformations, we obtain that if S is the ground- 
state configuration at (t, t + , t_, p e , p{), it is also the ground-state configuration at (— t, t + , t_, p e , p^. So 
the relative sign of n.n. and n.n.n. hopping amplitudes does not play any role in the model. 

Finally, let us note that Hamiltonian (fl|). in contrast to that where only n.n. hopping is present, is 
not invariant with respect to the joint hole-particle transformation for electrons and ions for any values 
of p e and pi. 
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Figure 2: The set of plaquettes (paths) over which the sums are taken in the ground-state energy 
expansion (j3J). 



3 The strong-coupling expansion of the ground-state energy 

Using the method of unitary-equivalent interactions ^7j, in the strong-coupling regime and at half-filling 
we can expand the ground-state energy Eg into a power series in t a t b + t c _. The result, with the expansion 
terms up to the fourth order (the fourth- order effective Hamiltonian) , i.e. a + b + c ^ 4, reads: 



E S (jl e ,fli) =4^ (Ve,^)+R (i) 
Es ] (fie, fit) = 



1 3 
- (/ii - fi e ) + -t 2 (t + + t. 



+ 



+ 



+ 



+ 



¥ - iV - iV (3t + + m+t - + 3e) + T 
it 2 , + -t A - -t 2 (2t 2 + + M-) - -t\ - -tit 2 

4 + 16 8 V + + ; 16 + 8 + " 

it 2 + -t 4 - -t 2 (2t 2 + t+tJ) - -t 4 - -t 2 ,t 2 

4" 16 8 V " + ; 16" 8 + ~ 



S xSy + 

J ( x ,y)z,+ 



~t 4 ~ -tH + t_ + — 1\ t 2 

8 8 + 16 + ' 



E3 22 X ^ ^ 2 2 X ^ 



16 



4,+ 



(x,y)4 



1 * ^ \ 1 a ^ v 3 o \ v 3 2 \ v 



+ 1 [t 4 + 2t%t„] s P2 + s P3 + ^V- J2 s p++ 

Pi Pz Pa 



16 

+ — t 2 t\ V S p + + —t 2 t 2 V , 

16 + L^i p 5 1 R - 



16 



(5) 



Here, the sign • in the (x, y)^. means a positive or a negative slope of an i-th order n.n. pair. The 
sets of plaquettes (or paths) P* are shown in Fig. |21 The sign • in the superscript of the path P' 
reflects the sign of the slope of the n.n.n. pairs in it. The remainder is independent of the chemical 

potentials and W, and collects all the terms proportional to t a t b + t c _, with a + b + c = 5, 6, The above 

expansion is absolutely convergent for sufficiently small t, t + and uniformly in A. In the special case 
of t + = £_ = t' and W = it was obtained in [THj . 

Let us note that the ground-state energy expansion (J3J) depends only on (^ — // e ). Hence, we denote 
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this difference as the unique chemical potential parameter /i. 

In Ref. [8], we have obtained the phase diagram of the isotropic model without n.n.n. hopping, 
i.e. for t+ = t- = 0. Here, our aim is to determine the influence of the n.n.n.-hopping anisotropy on 
the diagonal-striped phase Sdd (see Fig. |TJ). For this job, the value of the n.n.n.-hopping intensities, 
t+, t-, cannot be too large, in order to preserve the phase diagram up to 4th order. On the other 
hand, the n.n.n.-hopping intensities cannot be too small, in order to appear in the fourth-order effective 
Hamiltonian. In an attempt to satisfy the both requirements, we choose the smallest n.n.n.-hopping 
intensities t + , t_, i.e. such that they do not appear in the expansion terms of order smaller than four: 



_t 2 and £_ = a_t 2 . In this case, the effective Hamiltonian assumes the form: 



4 4) (/*) 



+ 



-fj, + -t (a+ + a 
2 4 v 



V 2 + ± t 4 
4 + 16 



X 

£_j S x$y + 



4 16 



* 4 + 



W 



S xSy + 



2,+ 



1 3 

4 1 G -+ I6 f 



^1 s x^y J r 



1 4 v 3 A ^ v 3 A ^ v 5 A ^ > 



5 

16 



p 2 



(6) 



i.e. the second requirement is satisfied. To answer the question put in the introduction, concerning the 
influence of anisotropy of n.n.n. hopping on the degeneracy of the phase Sdd, there is no need to consider 
the whole phase diagram. For t+ = t- = 0, we fix a point, well inside the domain of the diagonal-striped 
phase Sdd, sa y A* = an d W = —2t 2 + 9/2t 4 , i.e. to = 9/2 (see Fig. [TJ. Then, with the fixed point 
in (/i, H / )-plain, we introduce a n.n.n. hopping which does not change the ground-state configurations. 
Calculations show that Sdd has the minimal energy for a = \a + \ = |a_|, where —1/4 a 1/4 (we 
suppose that the difference between a + and a_ is not large, so they are of the same sign). Therefore, 
with our choice of n.n.n.-hopping intensities, the first of the above two requirements can also be satisfied. 
Eventually, we fix the values of n.n.n.-hopping intensities: a + = 1/8, a_ = 7a + , with 7 varying about 
1 (say, < 7 < 2). 

Since all the energy parameters, except the parameter 7 of n.n.n.-hopping anisotropy, have been 
fixed, the effective Hamiltonian © depends only on 7. In the following section, we examine how 
n.n.n.-hopping anisotropy, 7 7^ 1, influences the degeneracy of the diagonal-striped phase Sdd- 



4 Diagonal-striped phase versus n.n.n.-hopping anisotropy 

We use the m-potential method ^H] for constructing the phase diagram of effective Hamiltonian Eg (7) . 
For technical reasons, it is convenient to deal with such energies of configurations that are affine functions 
of the parameters of the effective Hamiltonian. However, the effective Hamiltonian ((HJ) contains the 
terms proportional to 7 and r y 2 . To get rid of nonlinearities, we replace 7 and 7 2 by two independent 
parameters d\ and di, respectively, with di, o?2 varying in the rectangle < d% < 2 and < c?2 < 4, in 
which the Hamiltonian is affine. After constructing the phase diagram in (di, G^-plane, we restrict it 
to the d 2 = d\ curve. 

To compare the energies of configurations, we rewrite Eg \d\,d2) as the sum, 

4 4) (^ 2 ) = !e f t\ (7) 
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over (3 x 3)-site blocks (called T-plaquettes) . The potential is of the form 



H. 



--(d 1 + i)s 5 + — j: 



3 

16 



S X Sy 



3 ^-v 3 5 -s—^n 



1/1 \ 1 

+ 32 V 16 2 + J ^ + 12 ^ SxSy + 
(x,y)2,- {x,y)3 



128^ ^ 128 ^ ^ 32 

p+ p~ P2 



(8) 



where s 5 is the central site of a T-plaquette. 

Unfortunately, the potential (jHJ is not an m-potential in the rectangle of considered values of d\ 
and d>2- Therefore, following |2lH |2*T| |2*2~] we introduce so-called zero-potentials. The zero potentials, 
satisfying 



(4) 
T 



0. 



(9) 



can be chosen in the form: 



(10) 



8=1 



where coefficients a, have to be determined in the process of constructing a phase diagram, and the 
potentials fe^, invariant with respect to the spatial symmetries of H , and fulfilling ©, read: 





—Si + Sg — 2S5, 






/Crp 


=S2 + s§ — 2s 5 , 








= s 3 + s 7 ~ 2S5, 








=S4 + S6 — 2s 5 , 








=Si^2 + SgSg — 


S4S5 


— S5S6 




= S2-S3 + S7S8 — 


S4S5 


— S5S6 




=SiS4 + SqSq — 




~ S 5 S 8 


/Crj-i 


=S3 S 6 + S 4-S7 — 


•S2S5 


~ S 5 S S 




=s 2 s 4 + s 6 s 8 - 




- S5S7 



Here, we label 1, . . . , 9, the sites of a T-plaquette, from left to right, starting in the bottom left corner 
and ending in the upper right one. In order to obtain the phase diagram, we have to compare the 
energies of all the possible T-plaquette configurations. The zero-potential coefficients a needed for this 
are given in Tab. ^ in the Appendix. We provide their values only at certain generating points, since 
we can assume that the coefficients a are affine functions of parameters (di,d 2 ). For more details on 
using the m-potential method see ^2] and references therein. 

The phase diagram of is shown in Fig. El The rectangle of considered points (^1,^2) breaks 
down into two domains. In the lower one, where 7 < 1 and £_ < t + , it is the phase S dd , with stripes 
being parallel to the direction of t_-hopping, that is stable. The analogous situation is in the upper 
domain, where 7 > 1 and £_ > t + : the stable phase, S dd , consists of stripes oriented along t + -hopping. 
At 7 = 1, we have the isotropic phase Sm whose configurations consist of S dd and S dd . 

So we see that, at least for the truncated effective Hamiltonian, switching on of a n.n.n. -hopping 
anisotropy reduces the rotational degeneracy of diagonal-striped phases: they become oriented in the 
direction of the weaker hopping. 
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o i d, 

Figure 3: The phase diagram of E^(di, (I2), with ^ di ^ 2 and ^ d 2 ^ 4. The rectangle breaks down 
into two regions: above d 2 = 1 the configurations S dd are the ground-state configurations, while S dd are 
replaced by S dd below d 2 = 1. The dashed line represents the condition d 2 = d\. The representative 
configurations (up to translations) of S dd , S dd are shown on the right. 

This result is similar to that described in ^2], where the influence of n.n. -hopping anisotropy on 
axial-striped phases was investigated. In that case, not only for a truncated effective Hamiltonian 
but also for the corresponding quantum one, it was proved that for any nonzero value of n.n. -hopping 
anisotropy the rotational degeneracy of axial-striped phases is reduced by making them oriented along 
the direction of the weaker n.n. hopping. Now in turn, the natural question is whether the conclusions 
we arrived at, concerning n. n.n. -hopping anisotropy, hold true for the quantum model, described by 
Hamiltonian (J3J). Applying the arguments presented in [2HI22], we can demonstrate (see for instance 
Ref. |23j). that the stable phases of the obtained above phase diagram remain stable for the model (jjj), 
but in some smaller domains. That is, if the remainder is taken into account, then there exist such 
a small to, that for t < t the phase diagram looks the same for the quantum model, excepting of some 
narrow regions (of width 0(t) in the scale of the fourth-order phase diagram shown in Fig. located 
along the phase-boundary lines. In our case that means that the breaking of the rotational symmetry 
occurs for 7 = 1 + 0(t), when the n.n.n.-hoping intensities are 0(t 2 ). Unfortunately, we cannot claim 
that any non-zero n.n.n.-hopping anisotropy reduces the rotational degeneracy of the quantum model, 
as it was the case for n.n. hopping (see [121 )• Here it seems, at least for small n.n.n.-hopping intensities, 
that there is certain critical value of (7 — 1|, above which the degeneracy of phase Sm is reduced. 

5 Summary 

We have considered the model of correlated spinless fermions, described by an extended Falicov- 
Kimball Hamiltonian. Quantum particles are allowed to hop between nearest-neighbor and next-nearest- 
neighbor sites. We have shown that a weak anisotropy of the next-nearest-neighbor hopping reduces 
the degeneracy of a diagonal-striped phase, it orients the stripes in the direction of the weaker next- 
nearest-neighbor hopping. 
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Appendix 

Here we present the zero-potential coefficients ol{ in the generating points of phase diagram shown in 
Fig. 01 

Table 1: Zero-potentials coefficients for the phase diagram shown in Fig. El 



Points (d\, d^) 


(0,0) 


(2,0) 


(0,1) 


(2,1) 


(0,4) 


(2,4) 


a\ 


83 
1536 


259 
3072 


79 
1536 


253 
3072 


67 
1536 


217 
3072 




25 
1536 


245 
3072 


29 
1536 


251 
3072 


41 

1536 


287 
3072 


a 3 


47 
1536 


331 
3072 


43 
1536 


325 
3072 


31 

1536 


289 
3072 




25 
1536 


245 
3072 


29 
1536 


251 
3072 


41 

1536 


287 
3072 


a 5 


53 
3072 


137 
8192 


19 
1024 


139 
8192 


3 
128 


91 
8192 


a 6 


53 
3072 


335 
24576 


19 
1024 


365 
24576 


3 

128 


91 
8192 


aj 


53 
3072 


401 
8192 


19 
1024 


403 
8192 


3 
128 


1117 
24576 


a 8 


53 
3072 


571 
24576 


19 
1024 


589 
24576 


3 
128 


139 
6144 


a 9 


19 
3072 


101 
2048 


19 
3072 


13 
256 


19 
3072 


1451 
24576 
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